Hyperacceleration in a stochastic Fermi-Ulam model 
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Fermi acceleration in a Fermi-Ulam model, consisting of an ensemble of particles bouncing between 
two, infinitely heavy, stochastically oscillating hard walls, is investigated. It is shown that the widely 
used approximation, neglecting the displacement of the walls (static wall approximation), leads to a 
systematic underestimation of particle acceleration. An improved approximative map is introduced, 
which takes into account the effect of the wall displacement, and in addition allows the analytical 
estimation of the long term behavior of the particle mean velocity as well as the corresponding 
probability distribution, in complete agreement with the numerical results of the exact dynamics. 
This effect accounting for the increased particle acceleration -Fermi hyperacceleration- is also present 
in higher dimensional systems, such as the driven Lorentz gas. 

PACS numbers: 05.45.-a,05.45.Ac,05.45.Pq 



In 1949 Fermi El proposed an acceleration mechanism 
of cosmic ray particles interacting with a time dependent 
magnetic field (for a review see |2(). Ever since, this has 
been a subject of intense study in a broad range of sys- 
tems in various areas of physics, including astrophysics 
H H plasma physics 0, atom optics H || and 
has even been used for the interpretation of experimen- 
tal results in atomic physics [lfj. Furthermore, when 
the mechanism is linked to higher dimensional time- 
dependent billiards, such as a time-dependent variant of 
the classic Lorentz Gas, it has profound implications on 
statistical and solid state physics E3- Several modifi- 
cations of the original model have been suggested, one 
of which is the well-known Fermi-Ulam model (FUM) 
\yX Eii, \\a\ which describes the bouncing of a ball be- 
tween an oscillating and a fixed wall. FUM and its 
variants have been the subject of extensive theoretical 
(see Re f. 1 131 and references therein) and experimental 
m eh Il7j studies as they are simple to conceive but 
hard to understand in that their behavior is quite com- 
plicated. A standard simplification [13j widely used in 
the literature, the static wall approximation (SWA), ig- 
nores the displacement of the moving wall but retains 
the time dependence in the momentum exchange between 
particle and wall at the instant of collision as if the wall 
were oscillating. The SWA speeds up time-consuming 
numerical simulations and allows semi-analytical treat- 
ments as well as a deeper understanding of the system 
[3 III H IH EH- However, as shown by Einstein in 
his treatment of the Brownian random walk j^, taking 
account of the full phase space trajectory (instead of the 
momentum component only) is essential for the correct 
description of diffusion processes. More recently, in the 
context of diffusion in the deterministic FUM, Lieberman 
et al have shown that one has to employ both canonical 



conjugate variables (position and momentum) in order to 
obtain the correct momentum distribution in the asymp- 
totic steady state • The present work shows that even 
in the absence of an asymptotic steady state the diffusion 
in velocity space is deeply affected by the location of the 
collision events in configuration space. 

The dynamical system in question consists of two har- 
monically driven infinitely heavy walls with an ensem- 
ble of particles bouncing between them. When a par- 
ticle collides with a certain wall, a random shift of the 
phase of the other wall, which is uniformly distributed in 
[0, 27r), occurs. The stochastic component in the oscilla- 
tion law of the wall simulates the influence of a thermal 
environment on wall motion and leads to Fermi acceler- 
ation El El El S El- It should be noted that 
although stochasticity can be introduced otherwise -for 
instance, via a random component in the angular fre- 
quency of oscillation- the random phase approach has 
become quite common as a method of randomization of 
the FUM and its modifications 

El El El, partially be- 
cause it is the only conceivable way to randomize the 
system without changing the energy of the moving wall. 

Despite the external randomization IS] the SWA does 
not provide an accurate description of the diffusion pro- 
cess; more specifically, the energy gain of the particle is 
substantially underestimated. For this reason we intro- 
duce in this Letter the so-called hopping wall approxima- 
tion (HWA), which takes into account the effect of the 
wall displacement. By means of this approximation it 
is made clear how the oscillation of the wall in configu- 
ration space affects the acceleration law of an ensemble 
of particles. Furthermore, the corresponding map allows 
analytical treatment and is as computationally efficient 
as the SWA and it enables us to calculate the evolution 
of the velocity distribution of the particles for long time 
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periods. 

The specific setup of the studied system is determined 
by the oscillation frequencies u>i and amplitudes Ai of the 
two walls (i = L, R) as well as the distance w between 
the walls at equilibrium. However, the dynamics does 
not depend on each of these parameters explicitly. It is 
therefore appropriate to introduce the relevant dimen- 
sionless quantities ej = r — Obviously, when the 
ratio x — meets the condition x > 1 (or x < 1) the 
original FUM is recovered. For the sake of simplicity the 
case x = 1 {Al = Ar = A) and r — 1 (ljl — = w ) is 
exclusively considered in the following. Using as a length 
unit the spacing between the walls w and as a time unit 
— the dynamical laws of the system can be derived in a 
dimcnsionlcss form: 



Vn 



-V n -i + 2u n ; u n = ecos(5t n + i n _i +i] n ) (1) 



where, U„ denotes particle velocity after the nth collision, 
u n for wall velocity on collision and rj n the random phase 
component. The time of free flight 5t n is obtained by 
solving the implicit equation 

X n -x + Vn-xStn = ±- +esm(St n + t n - X + r) n ) (2) 

Obviously, Eq. J3J) links the time of the nth collision to 
the position X n _\ of the particle in the previous colli- 
sion. The SWA simplifies the process on the basis of 
the assumption that the time of free flight, 5t n , depends 
only on particle velocity, i.e. St n = y 1 ■ If this ap- 
proximation is applied to the system of two oscillating 
walls, then it is possible to extract analytically the en- 
semble averaged velocity square (5V%) — (V* — V^_ x ) of 
the particle after integration over the random phase rj n : 

n 

(5V*) = 2e 2 . Given that <y„ 2 ) = £(^ 2 } + (V*), the 



i=l 



root mean square particle velocity is: 
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rms, SWA 



(n) 



2e 2 



(v 2 ) 



(3) 



Additionally, this quantity can be determined numeri- 
cally using the exact dynamical law given in eqs.fl},©. 
The calculations are performed on the basis of an ensem- 
ble of 10 4 trajectories with (V ) = e = j^. Cor- 
responding results are presented in Fig. 1 together with 
the analytical result ©, and show that there is a con- 
siderable difference between the acceleration rate of the 
root mean square (RMS) velocity given by the exact map 
and by the SWA. For a better understanding of how this 
difference originates, we improve the SWA by taking into 
account the impact of the displacement of the walls in- 
corporated in the exact dynamics. As particle velocity 
increases the time of free flight decreases, making it pos- 
sible to approximate the position of the scatterer at the 
instant of the nth collision with that at the (n — l)th 
collision. This approximation allows for an analytical 



evaluation of St n and defines the hopping wall approx- 
imation i.e. HWA. In this framework the time interval 
St n reads as follows: 
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where, 8t* r 



e[sin(t„_i+r)„)— sin(t„_2+?)n-i)] 



(4) 



is the correc- 



tion term to the time of free flight predicted by SWA. 

In order to derive (SV^) using eqs. Q), the following 
integrals have to be calculated: 



27T /■2-7T 



1 

i^ 2 



ecos 



(t n -i + St n + rjnj 



drindrin- 



(5) 

where j = 1,2. An exact analytical calculation of these 
integrals is not possible. However, for the set of parame- 
ters considered here, St n is much smaller compared to the 
other phase components. Therefore, we expand the r.h.s. 
of eq. © to the leading order of 6t n and then integrate 
over rj n and r\ n -\ which yields: 



e 2 cos 



2(V n -i) 



Therefore, we find: 



2e cos 



1 



e{V n - X ) 



2e z 



(6) 



In the limit of high particle velocities (V n -i) 3> 1, eq. © 
is simplified to (5V^) ~ 4e 2 , which is exactly two times 
the result obtained by neglecting wall displacement. Con- 
sequently the root mean square velocity as a function of 
the number of collisions is: 
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rms, HWA 



(n) = v /46 2 .n+ (V*) 



(7) 



The analytical result (JJJ based on the HWA is equally 
shown in Fig. 1. The above map can also be used to 
numerically simulate the acceleration process of the par- 
ticle, the corresponding results being presented in Fig. 1. 
In contrast to the static wall approximation which un- 
derestimates the acceleration of the particles, the HWA 
provides an accurate description of this process, indi- 
cating that the increased particle acceleration is due 
to the dynamically induced correlation between the po- 
sition and velocity of the oscillating wall on collision. 
This hyper-acceleration can be quantified by the ratio: 



Rh in) 



(SVJ) Exact 



2. However, it should be noted 



that the specific value of Rh depends on the character- 
istics of the oscillation law and more specifically on its 
turning points. For example, one can prove that for a 
piecewise linear oscillation law Rh is, in general, for any 
finite n different than 2 and only in the asymptotic limit 
becomes Rh (oc) = 2 251. 
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FIG. 1: Numerical results for Vrms of an ensemble of 10 4 
particles evolving in a FUM with two oscillating walls as a 
function of the number of collisions. Results were obtained 
by iterating the exact (circles) as well as the corresponding 
static (diagonal crosses) and hopping wall (upright crosses) 
approximative maps. It is noted that Vrma is measured in 
units of tow. Analytical results according to the SWA aver- 
aged over the random phase (solid line) as well as the analyt- 
ical prediction according to eq.JTJ (dash-dotted line) are also 
shown. 



The above analysis reveals the role of the fluctuations 
in the time of flight 6t n between successive collisions 
caused by the displacement of the scatterer. Despite 
the existence of an external stochastic component in the 
phase of the oscillating wall these fluctuations lead to a 
systematic increase of the acceleration of the particles. A 
simple explanation of the physical mechanism leading to 
the increased acceleration becomes possible by consider- 
ing the various configurations of the collision processes 
between the wall and the particles. Let us assume for 
a given velocity of the incident particle that the wall is 
moving in the same direction as the particle after passing 
the equilibrium position. The collision time due to wall 
displacement increases then, compared to the one assum- 
ing a wall fixed in space. In this case, the velocity of the 
harmonically oscillating wall is a decreasing function of 
time and therefore an increase of the collision time leads 
to a decrease of the wall velocity on the actual instant of 
the collision when compared to the static wall approxi- 
mation. This in turn leads to a lesser energy loss in the 
course of the collision. This reasoning holds equally in 
all other types of collision events, leading to the general 
picture of less energy loss or more energy gain when the 
wall displacement is taken into account. 

The focus of our attention now shifts to the probabil- 
ity distribution function (PDF) of the magnitude of the 
particle velocity and number of collisions n, p(\V\,n). It 
has been shown that the change of (V 2 ) as a function of 
the number of collisions n can be accurately described 



by a random walk in momentum space (V 2 ) oc n, pro- 
vided that the spatial motion of the walls is taken into 
account. Numerical as well as analytical treatments in 
setups similar to the present one suggest that p(|F|,n) 
is described by a spreading Gaussian 26, 27] . However, 
simulations with the exact map of eq.Q) yield the his- 
tograms in Figs. 2a, b /o(|V|, n) which show the PDFs cor- 
responding to snapshots for n — 5T0 4 , 5T0 5 collisions. It 
is shown that, even if the initial PDF is a Gaussian, with 
increasing time, it is transformed to a one-dimensional 
Maxwell-Boltzmann-like distribution since the set of ini- 
tial conditions leading to \ V\ <C 1 is vanishingly small for 
sufficiently long times. Using the hopping wall map the 
following analytical expression for the PDF is obtained, 
which describes the magnitude of the particle velocity for 
n > 1: 



p(\V\,n) = ±\V\e-l 



(8) 



where a = y 4e n +( v o ) _ j n pjg Sp 2a, b it is clearly seen 
that this analytical result predicted by the HWA accu- 
rately reproduces the exact behavior of the system. Fur- 
thermore, Fig. 2(c) shows the evolution of the mean value 
(\V\) as a function of the number of collisions n obtained 
numerically using the exact dynamics (open circles). For 
the sake of comparison we also show the corresponding 
analytical result: (|V|) = yf\o based on eq.lJSJ (solid 
line) . 
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FIG. 2: Numerically computed PDF for the magnitude of 
particle velocity using an ensemble of 10 4 trajectories and 
following the exact dynamics of eq.Q for: (a) n = 5 • 10 4 
and (b) n = 5 • 10 5 . In each case the analytical result given 
by eq.Q using the HWA is also shown (solid line). Finally, 
(c) shows the numerically obtained evolution of (|V|) as a 
function of n using the exact dynamics (open circles) as well 
as the analytical approximation based on eq.JHJ (solid line). 

The development of hyper-acceleration takes place in 
higher-dimensional scattering systems as well, such as 
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a time-dependent Lorentz gas consisting of harmonically 
oscillating circular hard scatterers on a triangular lattice. 
It is emphasized that in the time-dependent Lorentz gas 
system Fermi acceleration exists without any externally 
imposed randomization However, the absence or 

presence of a random component in the dynamics influ- 
ences the acceleration law. For example, if the oscillation 
axis is fixed and uniform throughout the lattice the ac- 
celeration law is V rms oc ti 1 / 4 . On the other hand, if the 
oscillation axis of the disks is randomly chosen on each 
collision, simulating the effect of thermal noise, the ac- 
celeration law becomes V rms oc n 1 / 2 , as in the ID FUM 
system |25j . In both cases, the static approximation un- 
derestimates the ensemble mean energy growth while the 
hopping approximation provides results much closer to 
those of the exact model. To illustrate this, in Fig. [3]nu- 
merical results for the random setup outlined above are 
presented. These are obtained utilizing the exact map 
|28| , as well as the corresponding hopping and static ap- 
proximative maps |22j. Consequently, it can be inferred 
that the development of hyper-acceleration is common to 
many driven dynamical systems, and features in any bil- 
liard which allows Fermi acceleration to develop. More- 
over, the understanding gained in the present investiga- 
tion helps open up the prospect of designing time-laws for 
the driving that provide a specified acceleration behav- 
ior for the ensemble of particles thus leading to a desired 
non-equilibrium i.e. time-evolving velocity distribution. 
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FIG. 3: Numerical results for K- ms of an ensemble of 10 4 
particles evolving in a triangular harmonically driven Lorentz 
gas with randomly chosen direction of oscillation on each col- 
lision, as a function of the number of collisions. Results were 
derived by iterating the exact maps (circles) as well as the 
corresponding static (diagonal crosses) and hopping (upright 
crosses) approximations. The parameters |A| = 0.01, U) = 1, 
w — 2.15, Vo — j-jg have been used in the numerical simu- 
lations, with A denoting the magnitude of the amplitude of 
oscillation, ui the angular frequency, w the spacing between 
the disk centers at equilibrium and Vo the initial particle ve- 
locity. Velocities are measured in units uiw. 
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